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Abstract. We study a generalization of the Kibble-Slepian (KS) expansion 
formula in 3 dimensions. The generalization is obtained by replacing the Her- 
mite polynomials by the q— Hermite ones. If such a replacement would lead 
to non-negativity for all allowed values of parameters and for all values of 
variables ranging over certain Cartesian product of compact intervals then we 
would deal with a generalization of the 3 dimensional Normal distribution. 
We show that this is not the case. We indicate some values of the parameters 
and some compact set in R 3 of positive measure, such that the values of the 
extension of KS formula are on this set negative. Nevertheless we indicate 
other applications of so generalized KS formula. Namely we use it to sum 
certain kernels built of the Al-Salam-Chihara polynomials for the cases that 
were not considered by other authors. One of such kernels sums up to the 
Askey- Wilson density disclosing its new, interesting properties. In particular 
we are able to obtain a generalization of the 2 dimensional Poisson-Mehler 
formula. As a corollary we indicate some new interesting properties of the 
Askey-Wilson polynomials with complex parameters. We also pose several 
open questions. 



1. Introduction 

In 1945 W.F. Kibble [2] and later, independently D. Slepian [3] have extended the 
Poisson-Mehler formula to higher dimensions, expanding ratio of the standardized 
multidimensional Gaussian density divided by the product of one dimensional mar- 
ginal densities in a multiple sum involving only the powers of constants (correlation 
coefficients) and the Hermite polynomials. The symmetry of this beautiful formula 
encourages further generalizations in the sense that the Hermite polynomials in the 
KS formula are substituted by their generalizations. 

Nice generalization of the Hermite polynomials emerged more than hundred years 
ago but only recently was intensively studied. The generalized Hermite polynomials 
are called the q— Hermite polynomials and constitute a one parameter family of 
orthogonal polynomials that for q = 1 is exactly equal to the family of the classical 
Hermite polynomials. 

We indicate that the function of 3 variables obtained with a help of so generalized 
KS formula has many properties of a 3 dimensional density. Let us call it fan- We 
know its marginals which are nonnegative. We are able to calculate all moments of a 
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supposed to exist random vector that would have this function as its joint density. In 
particular we can calculate the variance-covariance matrix of this " random vector" . 
This matrix is exactly the same as in the Gaussian case. Thus the question remains 
if function f^jj is really nonnegative for some values of parameter q (say — 1 < q < 1) 
and all values of correlation coefficients that make the variance-covariance matrix 
positive definite and almost all values of variables from certain product of compact 
intervals. It will turn out that it is not. We will indicate particular values of q, 
of correlation coefficients and a subset in R 3 of positive measure such that in 3 
dimensions so constructed generalization of the KS formula is negative on this set. 

Nevertheless we point out that it is worth to study described above sums since we 
obtain a nice and simple tool for examining properties of different kernels involving 
the so called Al-Salam-Chihara (ASC) polynomials. The kernels built of the ASC 
polynomials were studied by Askey, Rahman and Suslov in [5]. Such kernels have 
many applications particularly in connection with certain models of the so called 
q— oscillators considered in quantum physics. See e.g. j6], [7]. 

By studying the described above sums of the q— Hermite or ASC polynomials 
we obtain kernels that are different than those considered in [5]. Hence we obtain 
new results related to an important problem of summing kernels. 

One of these new results is a generalization of Poisson-Mehler expansion formula 
in the sense that the q— Hermite polynomials are replaced by the ASC ones. Of 
course the sum is different. Instead of the density of measure that makes ASC 
polynomials orthogonal (classical case) we get the density of measure that makes 
the Askey-Wilson polynomials orthogonal. We also analyze other non-symmetric 
kernels built of the ASC polynomials and sum them. As a by-product we point out 
in Remark [51 below the possibility of an interesting decomposition of the Askey- 
Wilson polynomials. 

The paper is organized as follows. In the next Section[2]we introduce all necessary 
auxiliary information concerning the so called q— series theory. In particular we 
introduce the q— Hermite and ASC polynomials and present their basic properties. 
In the following Section[3]we present our main results while less interesting or longer 
proofs are shifted to the last Section [5] We also include special Section [4] with open 
problems since not all questions that appeared when studying this beautiful object 
we were able to answer. 



2. Notation and auxiliary results 

Let us introduce notation traditionally used in the q— series theory, q is a param- 
eter. It can be real or complex, usually if complex, then \q\ < 1. We will assume 
however throughout the paper that —1 < q < 1. Having q we define [0} q = 0, [n] q 
= l + q + ... + q"-\ [n\\ = nlli [i] q , with [0] q \ = I, 

"1 =( F=W . . 

Mi q \ , otherwise 

It is useful to use the so called q— Pochhammer symbol for n > 1 : (a; q) n — 

1Yi=q (i" 01 ? 1 ) , with (a;q) = 1 , (oi, a 2 , . . • , a k ; q) n = U i=1 (a l ;q) n . 

Often (a; q) n as well as (eti, 02, ... , a^; q) will be abbreviated to (a) n and (a±, 02, . . . , a*) 

if it will not cause misunderstanding. 
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It is easy to notice that {q) n = (1 — q) n [n] \ and that 

{q) ", - , n > k > 
, otherwise 

Notice that for n > k > we have: [n] 1 = n, = n!, [™] = (™) (binomial 

cocfBcient), (a; l) n = (1 — a)" and [n] = 1 for n > 1, [ra] ! = 1, [£] = 1, (a; 0) n = 

f 1 if n = Ti , , r / s f 1 i/ i £ 4 
< . , ^ , • Let us also denote 1a(x) — { n ■ r a a ■ 

\ 1- a if n>l w \ if x £ A 

To define briefly and swiftly the one-dimensional distributions that later will be 
used to construct possible multidimensional generalizations of the Normal distribu- 
tions, let us define the following sets: 

S (?) = [-2/v/T - g, for |g| < 1 and S (1) = R. 

Further we define the following sets of polynomials: 
-the q— Hermite polynomials defined by the relationship: 

(2.1) H n+ i(x\q) = xH n (x\q) - [n]qH n -i(x\q), 

for n > with H-i(x\q) = 0, H (x\q) = 1, 

-the Al-Salam-Chihara (ASC) polynomials defined by the relationship: 

(2.2) P n+1 (x\y, p, q) = (x- pyq n )P n (x\y, p, q) - (1 - p 2 q n - 1 )[n] q P n -i(x\y, p, q), 

for n > with P_i (x\y, p, q) = 0, P (x\y, p, q) = 1, 

-the Chebyshev polynomials of the second kind defined by the relationship: 

U n +1 (x) = 2xU n (x) - Un-l (x) , 

for n > with U-x (x) = 0, U (x) = 1. 

Notice that H n (x|l) = H n (x) , H n (ar|Q) = J7 n (ac/2) and 

P„ (x||/, p, 1) = H n (-^M^j (1 - p 2 ) n/2 , where {i? n (x)} denote the 'probabilis- 
tic' Hermite polynomials i.e. monic polynomials that are orthogonal with respect 
to the measure with the density exp(— x 2 /2) / \/2tt. 
From Lemma [T] ii) , below it follows that: 

(2.3) P n (x\y, p, 0) = U n {x/2) - pyU n ^ {x/2) + p 2 U n ^ 2 (x/2) . 



The polynomials (|2.ip satisfy the following very useful identity originally formu- 
lated for the so called continuous q— Hermite polynomials h n (can be found in e.g. 
[19] Thm. 13.1.5) related to the polynomials H n by: 

(2.4) h n {x\q) = (1 - q) n/2 H n (7y=k) » « > 1. 

and here, below presented for the polynomials H n : for n, m > we have 

min(n,m) 

(2.5) H n (x\q)H m (x\q)= V W " [j] !i/„ +m _ 2j (x|g) . 

— ' 7 7 

One can also find in the literature (e.g. [16], [11]) the following useful formula: for 
m, n > we have 
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\j] q \H n -j (x\q) Hm-j (x\q) 



min(n.m) 

(2.6) H n+m (x)= ]T (-IYq® 

3=0 

that was originally formulated for the so called Rogers-Szego polynomials defined 
by: 



m 




n 


J _ 


9 


J. 



(2.7) 



x>, 



W n (x\q) = 

that are related to the continuous q— Hermite polynomials by: 

h n (x\q) = e in0 W n (e- m \q), 

for n > with x — cos 8. (See also [T9"].) 

To simplify notation let us introduce the following auxiliary polynomials of order 
at most 2 : 

(2.8) r k (y\q) = (l + q k ) 2 -(l-q)yq k , 

(2.9) v {x, y\p, q) = (l- p 2 ) 2 - (1 - + p 2 )xy + (1 - g)p 2 (ar 2 + y 2 ), 

(2.10) w fc (x, y\p, q) = v (x, y\pq k , q) , k > 1. 

It is known (see e.g. [9]) that the g— Hermite polynomials are monic and orthogonal 
with respect to the measure that has density given by: 



(2.11) 



/jV (x\q) 



{1)00 Vi - q 



271-^7-0 (x 2 \q) k=Q 
defined for \q\ < 1, x € M and 

(2.12) f N (x\l) 



1 



/2tt 



exp (~x 2 /2) 



Similarly it is known (e.g. from |12| ) that the ASC polynomials are monic and 
orthogonal with respect to the measures that for q E (—1, 1] have densities. These 
densities are given for \q\ < 1 by: 



(2.13a) fcN{x\y,p,q) 



2ny/r (x 2 \q) l} v k (x 7 y\p,q) 5(?) 



(,:) 



|p| < 1, x E E, y e S {q) and for q = 1 by: 

1 



few (x\y,p, 1) 



: exp 



(g - pvY 
'2(1 -p 2 ; 



It is known (see e.g. [H] formula 13.1.10) that for |g| < 1 

-n/2 



sup \H n {x\q)\ < W n (l\q)(l-qY 



(2.14) 



where W n is given by (|2.7|) . 

We will also need auxiliary polynomials {_B„ (a; |q) } Tl> 1 defined by the following 
3-term recurrence: 

(2.15) B n+1 (y\q) = -q n yB n (y\q) + q^ 1 [n] q B n . x (y\q) ; n > 0, 

with B_! (y|g) - 0, B (y\q) = 1. 
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In fact the polynomials {B n (y\q)} n >_i are equal to the polynomials \h n (y\q x ) } 
scaled and normalized in a certain way. The polynomials {h n (y|<7 -1 )} are known 
for a long time and were studied in [13] and [I]. However with the present scaling 
and normalization they were introduced in |12j where some of their basic properties 
were presented and their auxiliary role in finding connection coefficients between 
ASC and q-Hcrmitc polynomials was shown. One can show (see e.g. [T2]) that 
B n (x|l) = i n H n (ix) and 

( 1 if n = 0V2 
B n (x\0) = I -x if n = 1 . 
[ if n>2 

Further properties of these polynomials including their relationship to the q-Hermite 
polynomials are presented in |23j . 

Facts concerning the q-Hermite, ASC and polynomials B n , necessary for the 
derivation of the results of the paper, are collected in the following Lemma. 

Lemma 1. Assume that < q < 1, \p\ , \p^\ , \p 2 \ < 1, n, m > 0, x,y 1 z € S (q) , 
then 



i) 



ii) Vrt > I 



H n (x\q) = 



fc=0 



p n -«H n _ k (y\q)P k (x\y,p,q) 



P n [x\y,p,q) = 



fc=0 



p n - k B n _ k (y\q)H k {x\q) 



n r -i 

n 



E 

,;=o 



B n -j (x\q)Hj (x\q) = 0, 



m) 



^H n (x\q)H m (x\q)f N (x\q)dx = { ^, ^ ^ ™ , 



iv) 

S{q) 
V) 



Pn(x\y,p,q)P m (x\y,p,q)fcN(x\y,p,q)dx = ^ p2 ^ n] ^ ^ " f ™ 
-ffn /cat (x|y, p, g) = p n H n (y\q) , 



vi) 
vii) 



S(q) 



few {x\y, p t ,q) few (y\z, p 2 , q) dy = f CN (x\z, p 1 p 2 , q) . 



S(q) 



°° p n 
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viii) For \t\, \q\ < 1 : 

^ (ik) ^ _ j_ ^ wf{i\ q )r (t 2 ) c 

A 2 'Z^ 



OO 

§ («)< (?)< ~ w= 



convergence is absolute, where Wi (x\q) is defined by HI 
ix) 

Vx,y € S {q) :0 < C (y, p, q) < < — f 



f N {x\q) ~ ( p )l 

Proof, ii) was proved in [12] . i) is proved in [20] (formula 4.7) for the polynomials 
h n and Q n (x|o, &, g) related to the polynomials P n by the relationship Q n (x\a, b, q) 

= (1 — q) n ^ 2 P n ( ?f | ; 2a , \/&, <? ) ■ However one can also derive it easily from 



ii). iii) and iv) are known (see e.g. |19| ) for polynomials h n and Q n . Thus here 
they are presented after necessary adaptation to polynomials H n and P n . v) and 
vi) can be found in [9] , but their particular cases in different form were also shown 
in [8]. vii) is in fact the famous Poisson-Mehler formula which has many proofs. 
One of them is in [19] the other e.g. in [4] or [22]. viii) see Exercise 12.2(b) and 
12.2(c) of [19 . ix) was proved in [23J (Proposition 1 vii)). □ 

Let us remark, following [53], that 

fcN (y\x, Px,q) fcN (x\z, p 2 , q) 



(2.16) f A w(x\y,p 1 ,z,p 2 ,q) 



fcN (y\z,p 1 p 2 ,q) 



is the density of the measure that makes the re-scaled Askey- Wilson (AW) polyno- 
mials orthogonal. The AW polynomials mentioned in the paper are considered for 
certain complex valued parameters related to y,p 1 ,z,p 2 - For details see formula 
(2.5) of [33] . We will call the function defined by (|2"TT5]| the AW density. 

As mentioned in the Introduction the main object of this paper is a generalization 
of Kibble-Slepian formula. Traditional KS formula has a form of an expansion of 
the ratio of the non-degenerated n— dimensional Gaussian density divided by the 
product of its one dimensional marginals in the multiple series (in fact involving 
n(n — l)/2 fold sum) of the Hermite polynomials in one variable with coefficients 
that are powers of off-diagonal elements of variance- covariance matrix. 
We will analyze only the generalization of its three dimensional version. It is simple 
to express and general enough to expose interesting properties and applications. 

Namely we will analyze the following function given by: 

3 

(2-17) f3 D (x!, x 2 , x 3 \p 12 , p 13 , p 23 , q) = J|/jv (xi\q) 

i=l 

( 2 - 18 ) x E rfu P h P h , H »i O^g) H i+* (*2|g) H i+k (x 3 \q) , 

i,i,k>o L J ?' m «' [ V 

where |p 12 | , | Px3 1 > | P23 1 < Let us denote 

1 P12 P13 
P12 1 P23 i A = det P- 

Pl3 P23 1 
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We will assume that parameters p 12 , p 13 , p 23 are such that 

A = 1 + 2p 12 p 13 p 23 - p\ 2 - p\ 3 - p\ 3 > 0. 




and similarly for other pairs of indices (1,2) and (2,3) since A is obviously sym- 
metric in p 12 , p 13 , p 23 . 

The above mentioned assumptions concerning the parameters p 12 , p 13 , p 23 will 
be assumed throughout the remaining part of the paper. Similarly we will assume 
that ^ 6 S(?), i = 1,2,3 unless otherwise stated. Besides from Lemma [JJ viii) , 
ix) it follows that all considered in this paper series for \q\ < 1 are absolutely 
convergent, hence we will not repeat this statement unless it will be necessary. 

Remark 1. Let us immediately observe that when q = 1, that is when H n (x\q) is 
substituted by H n (x) and Jn (xi\q) by exp (— a;f/2) /V^n, then by Kibble-Slepian 
formula (see e.g. [3] Example 2)) for n — 3 , f%D (x\, x 2 , x^\p 12 , p 13l p 23l 1) is the 
density of the Gaussian distribution N (0, p). 

Besides let us notice that the function f 3 o has the same symmetry with respect to 
variables X\, x 2 , x 3 , p 12 , p 13 , p 23 as the density of the Normal distribution N (0, p) . 

Remark 2. In the literature functions somewhat similar to (|2. 17112. 18]) were consid- 
ered. In [10] (formula 4.13) sums of the form J2k>o J^\ Hn + k ( x \l) H k (y\q) Hk (z\q) 

for all ra's while in [15] (formula 2.5) sum of the form J2 m ,n,k ( q ) ( q ) Tqj W ™+k ( a \l) W 7 
( for W n see (]2.7I) ) were analyzed. Both these results refer to and are based on 
the beautiful results of Carlitz [T7j. The obtained sums are very complicated, ex- 
pressed through the basic hypergeometric function 6 05 or 3 2 . These results show 
that summing products of more than 2 g-Hermite polynomials is a very subtle 
and complicated task. On the other hand the results of these paper dealing also 
with sums of products of 3 g-Hermite polynomials follow subtle path in this area 
leading to results of relatively simple form. By the way our results also refer to 
above-mentioned results of Carlitz. 

Let us present some immediate remarks concerning the function /3D. 

Remark 3. i) 




hD (xi,X2,x 3 \p 12 ,p 13 ,p 23 ,q)dxidx 2 dx3 = 1, 



ii) 




fsD (xi,x 2 ,x 3 \p 12 ,p 13 ,p 23 ,q)dxidx 2 = f N (x 3 \q) , 




and similarly for other pairs (1,3) and (2,3). 



s 
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In particular we have: 



S(q)xS(q)xS(q) 



xihD {x\,x 2 , x 3 \p 12 , Pi 3 , /°23> l) dx x dx 2 dx 3 = 0, 



S(q)xS(q)xS(q) 



x\f 3 D (x 1 ,x 2 ,x 3 \p 12 ,p 13 ,p 23 ,q)dx 1 dx 2 dx 3 = 1 



and again similarly for the remaining indices 2 and 3. 
hi) 

/ Hd (xi,x 2 ,x 3 \p 12 ,p 13 ,p 23 ,q)dx 3 = fcN (xx\x 2 ,p 12 ,q) f N (x 2 \q) . 

JS{q) 

iv) If p 13 = p 23 = 0, then 

/3D {x 1 ,x 2 ,x 3 \p 12 ,p 13 ,p 23 ,q) = fcN {xi\x 2 ,p 12 ,q) f N (x 2 \q) f N (x 3 \q) 
and similarly for other pairs (1, 3) and (2, 3) . 

Proof. In all assertions we apply Lemma [T] hi) , the fact that H± (x\q) = x and 
H 2 (x\q) — x 2 — 1 and also formulae (|2.5[) and Lemma Q] v). iv) follows directly 
from Lemma [1] iiv) . □ 

As it follows from the above mentioned Remark that f 3 o is a serious candidate 
for the 3 dimensional density. It has nonnegative marginal densities equal to the 
densities /at and /n/cn- Moreover if f 3 o was a joint density of certain 3 dimen- 
sional random vector it would follow from Remark [3] ii) that variance-covariance 

[ 1 Pl2 Pl3 

matrix of this vector would be equal to p 12 1 p 23 

[ Pl3 P23 1 

To analyze its properties deeper we will need the following Lemma. 

Lemma 2. Let us denote j mJ , (x, y\p, q) = J^o ]$l H i+m (x\q) H l+k (y\q) . 
i) Then 



(2.20) 



7m,fe {x, y\p, q) = 7 n {x, y\p, q) Q m ,k (x, y\p, q) , 



where Q m .k * s o. polynomial in x and y of order at most m + k. 
Further denote C n (x, y\p x , p 2 , p 3 , q) = J2i=o [i\ q Pi~ l PlQ n—iii \X> y\P 3 ,Q) 
Then we have in particular: 
ii) 



]T(-i)^) 



k,m (V, x\p, q) = Q m ,k (x, y\p, q) 



p s H k _ s (y\q) P m+S (x\y, p, q)/(p 2 ) m + s , 



for all x, y g S (q) and q £ (—1, 1], 
Hi) 



Cn (x,y\p 1 ,p 2 ,p 3 ,q) = 



n r n 

n 



s=0 



H n _ s (y\q) P s (x\y, p 3 , q) p n x s p s 2 (p 1 p 3 /p 2 ) s / (pj) 
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Proof. Assertions i) and ii) are proved in |23j . Thus we will prove only iii). We 
have C n (x, y\p x , p 2 , p 3) q) = 

£?=0 [?] /l"Vi E;=0 ["71 ,// : V':;//, i , (»|?) (»|y, P3, ?) /(p§)<+J = 

£? =0 [:]H n - s (y\q)P s (x\y,Ps,q)pr s P S 2/(p 2 s) s EUo [;] q (-l) j <l®(PiPs/P2) j = 

Es=0 Sg-^n-s 0%iP3,9)Pl~ S P2 {PXPJP2) S I (Pl)s ■ 0n th<3 ^ W<3 haVe 

used the following identity (a) n = £™ =0 ["] q (~ l Y Q^a d - □ 

We get immediate observations: 
Corollary 1. For all n > 1 we /iawe: 



P„ (s|y, p, q) = ( P 2 ) n J2 n ^p'Hn-i (x\q) P, (y\x, p, q) I (p 2 ) i , 



ii) 

iii) 
iv) 
v) 



P n (x\y, p, q) fcN (y\x, p, q) dy = (p 2 ) H n (x\q) , 



S{q) 



C n (x,y\p 2 p 3 ,p 2 ,p 3 ,q) = p 2 H n (x\q) , 
C n (x,y\p 1 ,p 1 p 3 ,p 3 ,q) = p\H n (y\q) , 



C n (x,y\p 1 ,p 2 ,0,q) = ^2 



pr s p & 2Hn-s(y\q)H s {x\q), 



C n (x,y\0,p 2 ,p 3 ,q) = p 2 l P n (x\y,p 3 ,q) J (pl) n - 
C n (x,y\pi,0,p 3 ,q) = p'lP n {y\x,p 3 ,q) / (p 2 3 ) n . 



vi) 

C n (x,y\p 1 ,p 2 ,p 3 , 1) 



P1+P2- 2 PlP 2 P-j \ 

i-pI ) 



1/2 



x{p 2 - PiP 3 ) + y(Pi - P2P3) 
Vi 1 - pJKp! + p\ - 2P1P2P3) , 



Proof, i) see Corollary 2 of [23] . ii) We use previous assertion and Lemma[T]iv). iii) 
We have {p 2 pll P2) s = (pf) s , hence C n {x,y\p 2 p 3 , p 2 , p 3 ,q) = p% £" =0 ["] q H n - s (y\q) P s (x\y, p 3 , q) p 3 ~ 

f 1 if s = 

= p$H n {x\q) by Lemma [T]i). iv) (piP 3 /piP 3 ) s = < Q ^ s > q v ) First two 

statements are direct consequence of the assumptions that p 3 — or p 1 — 0. The 
third one follows the fact that p 2 (piP 3 / ' p-z) a = G°2 ~ Q t l PiP3) > which for p 2 

= is gives (— l) s q^ pfp 3 - Then we apply assertion i) of this Corollary. 
To get vi) first we notice that 



P n (x\y,p 3 ,l)/(l-p 2 3 y 



H„ 



x- p 3 y 



/ (1 - pI) 
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Then we apply the well known (see e.g. |18j ) formula for addition of Hcrmitc 
polynomials : 



V ( n ) a^Hn-i (0 Hi (9) = (a 2 + 6 2 )" /2 H n 



b 2 

with a = On, & = (P2 7 PlP , 3) . £ = y and 6* = a 7 p3j/ , . Next we observe that a 2 + b 2 = 

Pl+pj-2PlP 2 P3 an( J a £ _|_ foQ _ x (P2"PlP3) + VjPl-P2Pa) □ 
1— p 3 1 — p 3 

3. Main results 

Applying Lemma [5] to the function f 3 o we have the following Proposition: 
Proposition 1. i) 

(3.1) f 3D (x 1 ,x 2 ,x 3 \p 12 , p 13 , p 23 ,q) = fcN {x 3 \x!,p 13 ,q) f N (xi\q) f N (x 2 \q) 

x E TTl Hs ( X2 \l) Cs ( X l' X 3\Pl2,P23,Pl3,<l) > 

similarly for other pairs (1, 3) and (2, 3) , 
(3.2) 

f 3 D(xi, x 2 , x 3 \p 12 , p 13 , p 23 , q) = fcN {xi\x 3 ,p 13 ,q) fcN {x 3 \x 2 , p 23 ,q) In (x 2 \q) 
P1 2 (P13P23/P12) . 

] 9 ! (pi3 ; pi 3 ) s 

similarly for other pairs (1, 3) and (2, 3) . 

Proof. Lengthy proof is shifted to section [5] □ 

As an immediate corollary we have the following formula: 
Corollary 2. 

P12 (^13/^23/^12)5 



£ ^rS^ff Ps (X1 |X3 ' Pl3 ' 9) Ps ( * 2 '* 3 ' P23 ' g) 

/cw (xi \x 2 , p 12 , q) ^ p\ 3 {p 12 p 23 /p 13 ) h 



s=0 ^9 



/cjv(a;i|a;3,^i3,9) ^ [*],' (p? 3 > Pis) 



EP13 \P\2V23I P13^fc D \ t-> / 1 \ 

— — — 2 — 2 — P k (x 1 \x 2 ,p 12 ,q)P k (x 3 \x 2 ,p 23 ,q) . 



Proof. From Proposition [T] ii) we get 

P 12 (P13P23/P12 ) 
s . „ W g ! (Pi3.P3s) 



E TTT7 2 M ' f - ^ fta, 9) P. (*»!*», P 231 9) 
l s V \P131P23Js 

fcN {xi\x 2 , p 12 , q) fcN (x 2 \x 3 ,p 23 , q) f N (x 3 \q) 
fcN (xi\x 3 , p 13 , q) fcN (x 3 \x 2 ,p 23 , q) f N (x 2 \q) 

E %%Hrr^ fta. <?) ft (x 3 \x2, P23 , q) . 

fc>0 I- J ? \Pl3iP23ls 



Now we use the fact that: f C N (x 2 \x 3 , p 23 ,q) f N {x 3 \q) = fcN (x 3 \x 2 ,p 23 , q) f N (x 2 \q) . 

□ 
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Corollary 3. If p 12 = p 13 p 23 then 
i) 

Hd (xi,X 2 ,X 3 \p 13 p 23 ,p 13 ,p 23 ,q) = fcN (xi\x 3 , p 13 , q) f C N (#3 \x 2 , P23 , 0) /iV (x 2 \q) , 

ii) 

(3.3) 

^3 p r ^,_ „„ ^_ /cjv(Xl|X3,p 13 ,g) 

/cjv (a;i|a;2,Pi3p23:9) 



E 

s>0 L " J 9 

// p 12 = then 



ii,- 2 \ (^1 1^2,^13^23. ^ S (^3k2,P23.'?) = 

S U (PlsP^Js 



hi) 
(3.4) 



00 (— l) s gOp* 5 p s 

E r 1 1 ( 2 — irv 1 ^ (^1^3, Pi 3 , g) P* (a^ks, p 23 , o) 



l s \ q <- (Pl3>P2 3 ), 



E 



Pl3 



-P fe (xi|g)P fe (x 3 |a;2,p 2 3,g) 



/CJV (Xi|x 3 ,p 13 ,g) ^ W g !(P23)fe' 

Proof, i) We use Corollary[T]iii) and deduce that J2 s >o Wl Hs ( x ^\l) x a\pi 2 , P23> Pi3> <?) 

= J2 s >o ffi H » OM?) H s (x 3 \q) = fcN (x 3 \x 2 ,p 23 ,q) I f N (x 3 \q) (by Poisson-Mehler 

formula) or we notice that (Pi 3 P2 3 /Pi2) s when p 12 = P\ 3 p 23 is equal to for s > 1 
and use (|3.2p . 

ii) We use equivalent form of (|3.2[) that is 

Hd (xi,x 2 , x 3 \p 12 , p 13 , p 23 , q) = fcN {xi\x 2 , p X2 , q) fcN (x 2 \x 3 , p 23 , <?) /tv (x 3 \ q) 

P13 (P12P23/P1 



E 



-P s (xi|X2,Pi 3 P2 3 ,g)-P S (^3^2,P 23) 9) 



f^ Q Wg ! (Pl3»P2 3 ) s 

apply assumption, observing that then (Pi2P23/Pi3) s = (P23) s an d using the asser- 
tion i) of this Corollary. 

iii) We use the fact that p\ 2 (p 13 p 23 /p 12 ) s = ILLi ^12-9^^23) which for P12 
= equals to (— l) s q 1+ " +s ~ 1 pl 3 p 23 and the fact that P s (x\y, 0, q) — H s {x\q) . □ 

We have also the following remark concerning the ordinary, probabilistic Hermitc 
polynomials 

Remark 4. We have 



(3.5) 



/ (! ~ P13X 1 ~ P23) ~ (P12 ~ P13P23) 

V (i-p? 3 )(i-pI 3 ) 

(Pl2 - Pl3P23) S „ ( Xl~p 13 X 3 \ ( X 2 ~ p 23 X 3 \ 

^ S !(l-Pi 3 ) s/2 (1-P^ /2 "{V^pfs) S l«j 



exp 







1 


Pl2 


Pl3 


-1 


1 


P13P23 


Pl3 


-1 




Xl 




Xl x 2 x 3 ] 




P12 


1 


P23 




P13P23 


1 


P23 






X 2 


) 






. Pl3 


P23 


1 




Pl3 


P23 


1 






. X 3 





Proof. We use the fact that P, 



> n (x\y,p,l)/(p 2 \l) n = H n (-^=) /(I - P 2 )™ /2 - 



Then £~ , 



S=0 S !(l-pi 3 )./2(l-^3)./2"* I yiT^ 



,'!\ j ■'! -/-, j //s J \ can be summed by 
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the classical Poisson-Mehler formula yielding 

(3 ' 6) 7T=^ CXP( 20^3) } ' 

with r = (P12-P13P23) £ = ^-p 13 X 3 , = X 2 -p 23 X3 J t matter Q f e l emen _ 

\/(i-pf 3 )(i-pi3) V^ 3 ^ nA^pII 

tary calculus to check that the exponent of (|3.6p equals to the exponent of the right 
hand side of (l3"3i □ 



Remark 5. Notice that we get a new proof of the classical KS formula in 3 dimen- 
sions, since the left hand side of (|3.2[) is given by (|2 .17112 .15)) while the right hand 
side can be, in view of the Remark [4] written as 



V27r(l-p2 3 )(i_ p 2 3) 





1 


P13P23 


P13 


-1 


Xi 


Xi x 2 x 3 ] 


P13P23 


1 


P23 




X 2 




Pl3 


P23 


1 




X3 



( 1 ~Pl 3 )( 1 ~P23)~(Pl2-Pl3P23) 2 

(i-Pif3)(i-^ 3 ) 



x exp 







1 


Pl2 


Pl3 


-1 


1 


P13P23 


Pl3 


-1 




Xl 


Xi X 2 X 3 ] 




Pl2 


1 


P23 




P13P23 


1 


P23 






X2 






. Pl3 


P23 


1 




Pl3 


P23 


1 






. X 3 _ 





1 


Pl2 


Pl3 


-1 


Xl 




Xl x 2 x 3 ] 


P12 


1 


P23 




X2 


) 




. Pl3 


P23 


1 









27r (! - Pis)! 1 - P23) - (P12 - P13P23) 

which is the density of N (0, p) distribution. Thus (|3.2j) is really a generalization 
of 3 dimensional KS formula. 

We have also the following observation concerning the AW density. 

Remark 6. Following observation (|2.16p and Lemma [T]vii) we deduce that 

fcN (Xl \x 3 , p 13 , q) fcN (x 3 \x 2 , P 23 , <?) 



fcN (x 1 \x 2 ,p 13 p 23 ,q) 



fAw(x 3 \xi,pi 3l x 2 ,p 23 ,q)- 



Hence we have the following expansion of the AW density: 
(3.7) 

Pis 



f Aw{x 3 \xi, pi 3 ,x 2 , p 23 ) = Jcn (x 3 \x 2 ,p 23 ,q)^2 



. l*],' (Pfsffia). 



-P s (xi\x 2 , Pi 3 p 23 ,q)P s (x 3 \x 2 ,p 23 ,q) , 



which is an analogue of the Poisson-Mehler expansion formula interpreted as the 
one dimensional expansion given in Lemma [1] vii) . This result has been obtained 
by other methods in [24] , Besides from (|3.7j) it follows that: 



V?i > 1 



>S(q) 



(P„ (x 3 \x 2 ,p 23 ,q)- 



P13 fari f^ ( Xl \x 2 ,Pi 3 p 23 ,q))f A w (x 3 \xi,p 13 ,x 2 ,p 23 )dx 3 = 0. 
\pl3p23Jn 
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Thus there must exist functions Fi >n (x\, :E2|Pi3, P231 1) sucn that: 
V j > 1 : A 3 (aj 3 |xi,p 13 , x 2 , p 23 , q) = 

3 



^F hn {xx,X2\Px^p^q) (P n {x 3 \x 2 ,p 23 ,q) 



n=0 



Pl3 (P2z)n 



Pn (xi\x 2 ,p 13 p 23 ,q)), 



where {Aj}j >0 denote the Askey- Wilson polynomials with complex parameters 
related to x\, p 13 , x 2 , p 23 by the formulae (2.17)-(2.20) of [23] , 

As far as general properties of the function f 3 jj are concerned we have the fol- 
lowing formula that expresses function C n in terms of q— Hermite polynomials. 



Proposition 2. Vrt > 1 



[n/2] 



C n (x u x 3 \p 12 ,p 23 ,p 13 ,q) = , 2 \ V (-l) k q^ 

\Pl3Jn 



k=0 



n 




'2k 


2k_ 


Q 


k_ 



n-2k 

E 

i=0 



n - 2k 
i 



pi 3 ( 



n—i—2k 

P12 



; ( Pl3P23 gk 

V P12 



V P23 

Proof. Long proof is shifted to Section [5] 
Let us remark that when p 13 = P\ 2 p 23 then 



r;i 1 k k k I Pl2Pl3 Pl3P23 

[k\ q l -Pl2Pl3p23 . 

q V P23 Pl2 , 



Hi {x\\q) H n - 2k -i (x 3 \q) 



n—2k—i 



□ 



c n (xi,x 3 \p 12 , P23, P12P23, q) 

n—2k r 

E 



1 



[n/2] 



E (-i)**® 



a - 2k 
i 



(Pl2P2 3 )n fc=0 



n 




'2k 


2k_ 


9 


k 



Hi( Xl \q)H n . 



2k— i 



(x 3 \q) 



the formula obtained in [23] (3.2-3.3) in the context of Askey-Wilson polynomials. 

As stated in the Introduction the problem of non-negativity of the function f 3 £> 
for all allowed values of q and p's has negative solution. Above we indicated that if 
P12 = /°i3/ 9 23 ( or similarly for some other pair of indices) then f 3 o is positive for all 
— 1 < q < 1 and Xi,x 2 ,x 3 € S (q) . Now we will indicate some another relationship 
between p's and q so that for some values of Xi,x 2 ,x 3 £ S (q) the function f 3 o is 
negative. 

We have the following Theorem 

Theorem 1. Let \q\ < 1 and p 12 = qpi 3 p 23 , then there exists a set S C S (q) X 

S (q) x S (q) of positive Lebesgue measure such that for all (xi,x 2 ,x 3 ) G S, f 3 D is 
negative. 

Proof. Let us take p 12 = qp\ 3 p 23 and consider (|3.2p , then p\ 2 {p\ 3 p 23 / P\ 2 ) s = 
1 if s = 

-(!-?)pl3P23 if s = 1 ■ 

if s > 1 

Hence 

/3.o(a;i,a:2,:E3kPi3/>23> P13; P23> ?) = 

(! - 9)^13^23^1 ~ Pl3 X 3)(X2 - P 23 X 3 ) 



fcN (xi\x 3 ,p 13l q) fcN (x 3 \x 2 ,p 23 ,q) f N (x 2 \q) (1 



(1 - p 2 13 )(l - p 2 23 ) 
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Thus the sign of J 3 d is the same as the sign of 

1 _ (1 - g)Pl3P230l - Pl3 x 3 ){x 2 - p 23 X 3 ) 

(1-P 2 i3)(l-Pl 3 ) 

Equivalently f 3 o would be positive if for all Xi, X2, £3 € S (q) , \p 13 \ , |y0 23 1 < 1 and 

(l-pi 3 )(l-p 2 23 )>(l-q) 2 pl 3 p 2 23 
which comes from the positivity in (|2.19p we would have 

(! - P13X 1 - P23) > (! - 9)^13^23(^1 ~ Pis^sX^ - ^23^3)- 



Since = y/1 — qxi € [—2,2] the last inequality reduces to the following : 
(3-8) (1 - p\ 3 ){l - p\ 3 ) > P 13 P 23 (yi - Pi3V3)(y2 ~ P23V5) 

Let us select y\ = 2p 13 and y 2 = —2,023. The inequality now takes a form 

(l-P?3)(l-P2 3 )>Pl3P23(4-y 3 2 )- 

Now it suffices take p 2 3 , p 23 = 0.6, 1 > q > 1/3 and y 3 < 4 — |. On one hand one 
gets 

(l-p 2 3 )(l-/4) = .16> .6 2 

while on the other we have 

(1 - p 2 3 )(l - p 2 3 ) = (1 - .6) 2 =< .6 2 (4 - (4 - 1)) < /4p 2 3 (4 - y 2 ). 

Since the function / 3 £j is continuous in xi,x 2l x 3 hence there exists a neighborhood 
of points (2p 13 / (y/l — q, 2p 23 / y/l — q, 2y 3 /y/l — q) of positive Lebesgue measure 
such that this function is negative on this neighborhood. □ 

As a corollary we get the following fact. 
Corollary 4. 

, 39 n fcN{xi\x 3 ,p 13 ,q) (1- q)p 13 p 23 (xi- p 13 x 3 )(x 2 - p 23 x 3 ) 

fc'N(xi\x 2 ,qp 13 p 23 ,q) (1-Pi 3 )( 1 -P23) 

°° Pl3 (1 - 9V23) 

= ^ M ^<?P\ 3 P 2 23 ) S {1-P 2 23 ) PS (xlla;2 ' 9P13P23 ' 9) Ps (a;3l:E2 ^ 23 ' 9) ■ 



s =0 ^ 



Proof. Follows the (|3.2|) and the fact that for p 12 = qp\ 3 p 23 p\ 3 (Pi 2 P 23 / Pi 3 ) s — 
Ph Ms) a and further that = . □ 



Remark 7. Notice that we can take qp 13 p 23 instead P\ 3 p 23 in (|3.3[) getting fj^xi^q'pi^l 1) 

= E^=0 r.i t/Jll „2 \ (afl|*2,3Pi 3 P 2 3:9) P * ( x 3 |»2 , P 2 3 > «) ■ HenCe the POSitivity 

conditions for kernels of the form YlT=a as G°i> P2) ( x i\ x 2, Pi, q) Ps (x 3 \x 2 , p 2 ,q) 
are quite subtle. Besides comparison of (|3.3I) and (|3.9p can be the source of many 
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interesting kernels involving Al-Salam-Chihara polynomials. In particular for p l2 

= Q k Pl3P23 we have 

fCN(xi\x 3 ,p 13 ,q) A q k3 Pl3P23 (g 

' 2^ [ I I ? 2 N ( 2 \ P 3 \ X l F3, Pl 3 , 9) -Pi (^2 |X 3 , p 2 3> ?) 



Jew (a:i|z2,gVi3P23,9) £^ [jU (p? 3 ) 7 - (p| 3 ), 



EPl3 P23) g „ / , u \ -n r \ \ 

U\ Un2kJ n 2 \ ( n 2 \ P s [Xl\X2,q P 13 P23^) P s fa I *2 , P 23 , Q) ■ 
s=0 ^ V " Pl3P23>s \P23)s 

4. Open Problems 

(1) Consider the KS formulae for higher (than 3) dimensions. Can it be non- 
negative for all allowed values of the parameters p and all values of the 
variables form S (q)? 

(2) Since for q = both densities fcN and /at and ASC polynomials (compare 
(|2.3[) ) are very simple one can hope to obtain exact formula for function 
./3.D- One can deduce that then f 3 £> divided by the product of the one 
dimensional marginals should have the form of the ratio of quadratic forms 
in 3 variables and the product of three quadratic functions of every pair 
of variables. Thus one can hope to find a set O in ( — 1, l) 3 such that if 
(P121 Pi3; P23) S O then f 3 n is a 3— dimensional density. 

(3) What about the cases p 12 = q k P\ 3 p 23 for fc > 1. Do we get non-positivity 
for all fc > 1 or are there some fc's for which it is nonnegative? 

(4) Numerical simulations suggest that for say p 12 = one can find q, P\ 3 , P23) 
x±,X2,X3 € S (q) such that one of the kernels (and consequently both) 
given in (|3.4|) is negative. What would be the proof of this fact? More 
precisely what would be range of parameters q, p 13 , p 2 3 an d subset A C 
S (g) 3 of positive Lebesgue measure such that for x±,x 2 ,x 3 G A the kernel 
is negative. 

(5) What about cases \p 12 \ > p 13 p 23 ? 

5. Proofs 

Proof of Proposition [7J First of all let us notice that due to inequality (|2.14l) and 
assertion viii) of Lemma [1] series defining function f 3 £> is absolutely convergent. 
Now to prove i) we apply Lemma [2] i) to (I2.17P and (|2.18[) and remembering that 
7o : o v\P> q) = fcN (x\y, p, q) /f N (x\q) we get: 

/ 3 d {xi,x 2 ,x 3 \p 12 ,p 13 ,p 23 ,q) 

P J i2P k 23 „ _,„^srti3zT {xilq)Hi+k (x 3 \q) 



i=l j,k>0 Ui 1' 1 V i>0 [ V 



= fcN (xi\x 3 , p 13 , q) f N (x 3 \q) f N (x 2 \q) ^ r Ttf 1 Hj+k ^ 2 ® j k ' x ^ Pl3 ' ^ ' 

j.k>o [Jl r [K W 

Now changing order of summation and substituting j + fc — > s, j— > fc we get 
f 3D (x!,x 2 ,x 3 \p 12 ,p 13 , p 23 , q) = fcN (x 3 \xi,p 13 , q) f N (xi\q) f N (x 2 \q) 

x E TTl Hs ( X2 \^ Cs {xi,x 3 \p 12 ,p 23 ,p 13 ,q) ■ 



1(5 



PAWEL J. SZABLOWSKI 



ii) We apply ([2^]) , Then 

3 

f3D(xi,x 2 ,x 3 \p 12 ,p 13 ,p 23 ,q) = Y[fN (xi\q) T~] 

i=l i,?",fc>0 ' ' 



P l2P23Pl3 

u*>o W 



H i+3 {xi\q)H j+k (x 2 \q) 



E 

n>0 



i 




k' 


n 


q 


n 



} q \ (-1)" g(5)fli_„ (aj 3 k)^fc-n (as 3 |«) 



00 00 



„k—n Jt—ri 



E E - n] 3 ! [fc 1 ! n l ^ 3 ' 9 ^ i/fe "" ^l 9 - 1 jffj ' +fe ^ 2 ' 9 ' ) 



_n fc=Ti 



3 00 n n 00 j 

1=1 n=0 



00 00 



E E [ffffr^ {x3lq) Hi+n+ i {xilq) Hk {X3lq) Hj+n+k {X2lq) ■ 



i=0 k=Q L '1' 1 >q 



Now we use quantities: j m k (x, y\p, q) = Y,iLa ^H i+m (x\q) H i+k (y\q) and apply 
(|2.20p and Lemma [2] iii) . We get then 

f3D(xi,x 2 ,x 3 \p 12 , p 13 , p 23 ,q) = JJ/jy(xilg)2j(-l) n gw 23 " 



1=1 



n=0 



x E r7n 7 o,™ +J (»3,a;i|pi3, g) 7o,«+j (^3,^2|p 2 3> g) 
= jew (xi\x 3 , p 13 ,q) Jcn (x 2 \x 3 , p 2 3><l) In (x 3 \q) 



V ( l) n a ( " ] p23Pl3 V ^ 

9' j=o u, i 



n=0 



, ^r-i - , /51 P «+J ( x i\ x 3,Pi3,q)Pn+j (x 2 \x 3 ,p 23 ,q) 

n '" ,=n Wl 3 Jn+j IP237,i+j 



E 



/CAT (Xi I S3 , p 13 , g) Jew (z 2 I X 3 , p 23 , g) / N (X 3 I g) 
1 



-P s (xi\x 3 , p 13 , q) Ps {x 2 \x 3 , p 23 , g) 



E 

n=0 



„_r W 9 !(P 2 3 ) S (^ 3 ) S 
(-1)9(2) (p 13 p 2 3)"Pi2™ = /cat {xi\x 3 ,p 13 ,q) fcN (x 2 \x 3 ,p 23 ,q) f N (x 3 \q) 



E 



P12 (^13^23/^12). 

jM,!(P?s).(4.), 



P s (a;i|a;3,p 1 3,g)P s (x 2 |z 3 , p 23 , g) 



□ 



Proo/ of Proposition® We start with C n (x, y\p 1 ,p 2 , p 3 ,q) = J27=o ["] q Pi 'piQn-hi (x, y\p 3< q) 
= ELo [H\ Hn -* (y\l) p s ( x \y>P3,<l) Pi~Vl {PiP3/Pi) s / (pI) s and a PP!y LemmaQ] 
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ii). 

n 

C n (x 1: x 3 \p 12 ,p 23 ,p 13 ,q) = E 



s=0 



p^ s H n _ s (x 3 \q)pt 3 

q V P23 



P s {xx\x 3 ,p 13 ,q) / (p\ 3 ) t 



1 - 

(Pia)n s =q 

1 n 



p™ 2 s H n _ s (x 3 \q) p s 23 



(Pl3)n s=0 



P12 S Hn-s (x 3 \q) p s 23 



P12P13 \ 
P23 /, 

P12P13 



P s (xi\x 3 ,p 13 ,q) (p%q s ) n _ s 



pU s ) 



1 n r "i n 



n — i 
s — i 



P23 

p^ s H n - s {x 3 \q) p s 23 



E 



Pi 3 l Bs-z (x 3 \q)Hi (xi\q) 



P12P13 \ 
P23 A 



Pfa* B '-i (x 3 \q) 



1 - 



-Hi Oik) x 



E 

j=0 



n — i 
3 . 



p^H^ (x 3 \q) (%> (Pll£±l) (p\ 3 q^) p\ 3 B, (x 3 \q) 



V P23 J 



i+3 



And further using formula 

L(n+m)/2J 

H m {x\q)B n (x\q) = (-l) n q^) E 

i=0 

proved in |23j (Lemma 2 i)) we get 



n 




n - 


- m — i 


i 






i 











1 

C n (xi,x 3 \p 12 ,p 23 ,p 13 ,q) = 2 E 

\Pl3Jn 



n r -1 

n 



i=0 



[f] !g-*(»-*)ff wfm _ ai ( a: |g) 



E 



n — « 
. 3 . , 
L(«-0/2J 

x E 

fe=0 



n— 2- 

P12 



V P23 /i+j 



n — i — k 
k 



[k] q \q- k ^H lMk {x 3 \q). 



Now keeping in mind that 



for j < k we split hrst internal sum into two 



sums : one form to [_(n — i)/2\ and second from \_(n — z) / 2J + 1 to (n — i). 

L(n-i)/2J r .., 

W(^) (Pis^l^^VisC-iF? 

V P23 /i+7 



*E " E 



2 

13/n j = Q L " J q 



3=0 



n — 1 
3 



n—i— 7 

P12 



E 

k=0 



f 




n 


— i — k 


A 


Q 




k 



(n-i) 

E 

i=L(n-0/2j+l 



n — i 
3 . 



n n-i-j J+j j P12P13 \ 
Pl2 P23 1 I 

7 \ P23 / t 



[k] Q \q- k ^ H lMk (x 3 \q) 

(pl 3 q^) n _^p{ 3 (-iyq& 



+3 



E 

t— 



f 




n 


— i 


k_ 


Q 




k 



[ k } q }q- k 0-k) Hn _._ 2k {x3 \ q)) . 
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We have further after changing the order of summation. 



1 " 

C„ (xi,x 3 \p 12 ,p 23 ,p 13 ,q) = 2 



L(n-i)/2j 

x( E 



[n-i] q \ 



n—i—2k 



Hi (x!\q) 



(x 3 \q) 



L(n-i)/2j 



E ?- fe(j - fe v 13 (-i) j 5^ I — ^ 



„' [n-t- j],! V P23 



L(n-i)/2j 

+ E 



[n-i\ q \ 



(n-i) 

x E 

j=Y(n-i)/2\+l 



[k] Q \[n-i-2k] q \ 
I 



H 



n—i—2k 



(a* k) 



[j -k] q \[n-i- j] q \ 



(ph) 



n r ~\ 

n 



Hi ixi\q) 



L(n-i)/2j 



'( E 



[n-i] 9 ! 



,„ [*],![n-i-2*],! 



n—i—2k 



(x 3 \q) 



(n-i) 

><E 

j=k 

Hence 



n — i — k 
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Now we use formula £" =0 (-if qd> [™] q {a) i V (abq l ) n i = (b) n proved in ([23]) 
(Lemma 1 ii)) with a — Pl p 2 ^ 13 q l+k , b = p 13 p 23 / 1 p\i getting 



1 
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Finally change the order of summation and using on the way an obvious property 
01 ( a ) n = ( a ) j ( a< ? : ')„_ for ever y < j < n 
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